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We present a theoretical study of Bose-Einstein condensation in highly anisotropic harmonic 
traps. The bosons are considered to be moving in an optical lattice in an overall anisotropic harmonic 
confining potential. We find that two-step condensation occurs for lattice bosons at much reduced 
harmonic potential anisotropy when compared to the case of an ideal bose gas in an anisotropic 
harmonic confinement. We also show that when the bosons are in an isotropic harmonic confinement 
CM but with highly anisotropic hopping in the optical lattice, two-step condensation does not occur. 

We interpret some of our results using single boson density of energy states corresponding to the 
1} , potentials faced by the bosons. 
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I. INTRODUCTION 



The experimental observation 1-3 of Bose-Einstein (BE) condensation in confined bose atom clouds have launched 
extensive experimental and theoretical studies of this phenomenon and the various properties of the condensates of 
free bosons 4 and lattice bosons 5 " 10 . In three dimensional non-interacting free 7 or lattice bose systems 11 in isotropic 
harmonic traps, the BE condensation is accompanied by a peak in the specific heat at the condensation temperature. 
However, not long after the initial discovery of bose condensation in harmonic traps, Druten and Ketterle (DK) 12 
found that bosons in highly anisotropic three dimensional (3D) harmonic traps show a qualitatively different behavior 
of the specific heat. Their theoretical calculations showed that the peak in the specific heat of this system is not 
i-^j , at the bose condensation temperature, but in a higher temperature range in which the bosons, as the temperature 
decreases, are progressively transferred from the tightly confined dimensions to loosely confined dimension. Thus, in 
this system, the bose condensation occurs in two steps: first the aforementioned transfer and second the condensation 
, ^ , ■ of all bosons in to the overall ground state. The calculations of DK are for free bosons in anisotropic 3D harmonic 
traps. The earlier 13 and the later works 14 are also for free bosons with anisotropic box and anisotropic 3D harmonic 
(f) i confinements, respectively. 

The two-step condensation is a phenomenon special in at least two aspects. The first aspect is that the condensation 
of the bosons into the overall ground state occurs without any signature in a thermodynamic property (the specific 
heat). The second aspect is that it involves a dimensionality cross-over in a higher temperature range. Both of these 
aspects are of fundamental importance, and hence an investigation of this phenomenon for lattice bosons should be 
of significant interest. Now, when an additional optical lattice potential is applied to bosons in a harmonic trap, it is 
(N ■ known that significant changes occur in the single boson energy density of states 11,15,16 (DOS). How these changes alter 
the two-step condensation phenomenon is an issue of fundamental importance and contemporary relevance. Further, 
considering that there is a dimensionality cross-over involved, it should of interest to know under what conditions a 
lattice bosons system in an anisotropic harmonic trap is truly one-dimensional as far as thermodynamic properties are 
concerned. Motivated by these considerations, in this paper we present a theoretical study of the condensation of lattice 
bosons in anisotropic harmonic traps. Among other results, we show that a two-step condensation occurs for lattice 
bosons at anisotropies much smaller than that required for free bosons. This finding may lead to an experimental 
study of this phenomenon in optical lattices. Further, since it is by now relatively easy to experimentally change 
the boson hopping in different spatial directions in an optical lattice, we also study a case of lattice bosons with 
anisotropic hopping in an isotropic harmonic confining potential. One may argue that this is equivalent to bosons 
with isotropic hopping in an anisotropic trap and hence two-step condensation would occur in this case as well. Our 
calculations show that two-step condensation does not occur in this case. Our focus in this paper is the effects of 
the optical lattice potential on the two-step condensation, not the effect of boson-boson interactions. The strength of 
boson interactions is controllable in optical lattices 9,17 , and hence one can experimentally reach the ideal gas limit to 
wheih our results are applicable. 

This paper is organized as follows. In the next section, we revisit the system of free bosons in an anisotropic 
harmonic trap. We complement the analysis of DK by new features we find in the single boson density of energy 
states (DOS). In Sec. II, we analyze two-step condensation of lattice bosons in anisotropic harmonic traps. In Sec. 
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Ill, we study the bose condensation of lattice bosons with anisotropic hopping in an isotropic harmonic trap. The 
conclusions are given in Sec. IV. 



II. TWO-STEP CONDENSATION OF FREE AND LATTICE BOSONS 



A. Free bosons in an anisotropic harmonic potential 



In this section we consider boson atoms in an anisotropic harmonic confining potential {K x x 2 + K y y 2 +K z z 2 ). The 
single boson energy levels are E(n Xl n yi n z ) ~ (n x + l/2)fku x + (n y + l/2)huj y + (n z + l/2)huj z , where n x ,n y ,n z = 0, 1, 
2.. .co, 2irfi is the Plank's constant, and 0J x , y , z = \J2K^~ y ~ z Jm in which m is the atomic mass. We use K x <sc.K y = K z 
so that the potentials, from the center of the trap, along the y and z directions are much steeper compared to that 
along the x direction. In order to exhibit the two-step condensation in this system, we calculate the ground state 
occupancy (No), the occupancy in the loosely confined direction (Nuj), and the specific heat by determining first 
the chemical potential (fi) from the bosons number equation N — _ _ X^=o N(E(n x , n y , n z )), where 
N(E) = l/(exp[/3(E - fi)} - 1), = l/k B T, T the temperature, and then calculating N = N(E(0,0,0)) and 
Nir> — X^=o N{E(n X} 0, 0))). The specific heat is obtained from the temperature derivative 11 of the total energy 
(Etot = X^= X^=o Z)£°=o N(E(n x , n y , n z )) E(n x ,n y , n z )). We evaluated N , N 1D , and C v numerically to obtain 
exact results. The results of our calculations for 10000 bosons in harmonic traps of varying anisotropics are shown in 
Fig. 1. The condensation temperature 11,18 (To), mentioned in this and other figures, is obtained by setting iV = 
and \x = £7(0,0,0) in the bosons number equation. For low anisotropy, the peak in C v is clearly associated with the 
growth of the bose condensate fraction (Nq /N) . With increasing anisotropy, the peak in C v becomes associated with 
the dimensionality cross-over in which bosons are transferred from the tightly confined directions (y and z) into the 
loosely confined direction(x). To complement the analysis of DK, we calculated the single boson energy DOS for these 
systems. First we recall that the analytical approximation to the DOS of a boson in a D-dimcnsional anisotropic 
harmonic trap 7 is 



From exact numerical calculations, we find that the DOS for highly anisotropic traps, shown in Fig. 2a, show features 
qualitatively different from the analytical approximation. The DOS is found to have one dimensional (ID) character 
(flat regions). In the very low temperature range, only the lowest flat region is occupied by bosons so that the C v 
shows ID character as shown in the inset of Fig. 1. This is consistent with the higher temperature cross-over in which 
the bosons are transferred from the tightly confined directions to the loosely confined direction. We also note that 
ID features are absent in the DOS when the anisotropy is comparatively smaller as shown in Fig. 2b. In the next 
section, we consider consider two-step condensation of lattice bosons. 



In this section we consider boson atoms in a three dimensional simple cubic (sc) lattice in an overall anisotropic 
harmonic confining potential. The system Hamiltonian is 



where t is the kinetic energy gain when a boson hop from site i to its nearest neighbor site j in the optical lattice, 
c\ is a boson creation operator, V(i) = (K x x 2 + K y yf + K z zf) is the potential at site i, n>i = c\ci the boson number 
operator, and fi the chemical potential. We numerically diagonalize this hamiltonian to obtain the energy levels (E{) 
of a lattice boson for different values of the anisotropy of the harmonic trap with k x <^ k y — k z , where k a = K a a 2 , 
a = x,y, z and a is the lattice constant of the sc optical lattice. The chemical potential and the boson populations in 
different energy levels are calculated at a temperature T, using the boson number equation: N = J27Lo N(Ei), where 
Eo and E m are the lowest and the highest single boson energy levels and N(Ei) = l/[exp[(E[ — /x)/fesT] — 1]. The 
specific heat is calculated from the temperature derivative of total energy (E tot = Y^iLn N(Ei)Ei). The lattice size, 
we use, is large enough that the results are free of finite size effects. The results for the temperature dependence of 
the bose condensate fraction (N /N), fractional boson population in the x-direction (Nm/N), and the specific heat 
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B. Lattice bosons in an anisotropic harmonic potential 




(2) 
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for lattice bosons in harmonic traps with different anisotropics are shown in Fig. 3. The two-step condensation is 
clearly seen in highly anisotropic traps (dash-dot and dotted lines), while it is absent for small anisotropics. Further, 
on comparing Figs. 1 and 3, we can immediately infer that the the two-step condensation of lattice bosons is possible 
for much smaller anisotropies compared to the free bosons in anisotropic harmonic traps. We also find that most of 
the lattice bosons are confined in the x-direction at much higher temperatures compared to the case of free bosons 
even though the anisotropy is much smaller for the former compared to the latter. This is clear if we compare Nid/N 
in both cases, as shown in Figs. 1 and 3, for T > 2.5Tb- In order to get some insight into the two-step condensation 
occurring for smaller anisotropies for lattice bosons, we calculated the single boson energy DOS for these systems. 
The results of our calculations are shown in Fig. 4. With increasing anisotropy there is a dramatic change in the 
single particle DOS (Fig. 4). Except for very low energies, the DOS for low anisotropy is very similar to that for the 
isotropic case. With increasing anisotropy, the DOS develops features similar to that expected in ID. It has peaks 
separated by broad regions of comparatively very small values of DOS. Since the DOS shows strong ID features for 
small anisotropies compared to the free boson case (without the lattice potential), the ID behavior in Nid/N sets in 
for comparatively smaller anisotropies of the harmonic traps. 

C. Dimensionality cross-overs of lattice bosons 

In this section, we investigate the dimensionality cross-over in some more detail. The specific heat, being a very 
sensitive function of dimensionality changes, is an useful property to study for this purpose. As mentioned earlier, for 
bosons in a highly anisotropic harmonic trap, the peak in the specific heat occurs in a temperature range in which the 
bosons are transferred from the tightly confined directions to the loosely confined direction (here x-direction passing 
through the center of the trap). With lowering of temperature, when the transfer to the single ID chain along the 
x-direction is complete, the system would show ID behavior. To visualize this cross-over, we have plotted the specific 
heats of lattice bosons in 3D anisotropic traps and in a ID trap in Fig. 5. It is seen that at low temperatures the 
specific heats of the anisotropic systems completely coincide with that of the ID chain showing that the anisotropic 
3D systems have become purely ID. The departure from the ID behavior starts at a temperature which is higher for 
larger anisotropy (Fig. 5). However, this temperature is much lower than the temperature where the specific heat 
shows the peak and is also lower than the condensation temperature (To) for anisotropy as high as k y /k x = 2 x 10 5 
for lattice bosons (see Fig. 5). It is also found that above this temperature (i.e., where the departure from the ID 
behavior starts in C„), a very small change in N\d produces a very large change in C v . This is also clear from Fig.l 
(for free bosons) and from Figs. 3 and 5 (for lattice bosons) on comparing the temperature dependences of C v and 
Nid in the above mentioned region. 

III. LATTICE BOSONS WITH ANISOTROPIC HOPPING IN AN ISOTROPIC HARMONIC 

POTENTIAL 

In this section, we investigate if the properties of lattice bosons with anisotropic hopping in an isotropic harmonic 
trap is similar to those with isotropic hopping in an anisotropic harmonic trap. We have studied the effects of 
anisotropic hopping on C v and condensate fraction in isotropic harmonic traps. The spatial direction dependent 
hopping we use are such that t x 3> t y = t z . The results of our calculations for various anisotropy ratios are shown in 
Fig. 6. We do not find a two-step behavior in this system. The specific heat is suppressed with increasing hopping 
anisotropy. The system shows anisotropic 3D behavior in all the properties studied even for large hopping anisotropies. 
The reason is that even though the hopping is anisotropic, since the trap is isotropic, one will have a large collection 
of ID chains, symmetric about the central one, coupled by the small inter-chain hopping. At the bose condensation 
temperature, the bosons condense in to the central region of the trap. There is no higher temperature cross-over 
region as in the case of lattice bosons with isotropic hopping in a highly anisotropic harmonic trap. 

IV. CONCLUSIONS 

In this paper we presented a theoretical study of two-step condensation of lattice bosons in anisotropic harmonic 
confining potentials. We find that, even though significant changes occur in the single bosons DOS, two-step conden- 
sation occurs for lattice bosons and it occurs for much smaller anisotropies of the harmonic potential compared to 
that for free bosons. This may lead to experimental observation of these type of condensations. We find unusual ID 
features in the one boson DOS both for free and lattice bosons in highly anisotropic harmonic traps. These DOS's 
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are equally well applicable to fermions. Using C v as the indicator, we also investigated the dimensionality cross-overs 
in the system. We find that the system shows purely ID behavior in a much lower temperature range compared to 
the temperature range in which the peak in C v occurs. Further, we investigated if two-step condensation occurs for 
lattice bosons with anisotropic hopping in an isotropic harmonic trap. We find that a two-step condensation does not 
occur in this system. 
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FIG. 1. The variation of the specific heat (C v ), bose condensate fraction (No/N), and the fractional population in the 
a;-direction (Nm/N) with temperature for free bosons in anisotropic harmonic traps. To is the condensation temperature as 
defined in the text. Solid line: K x = 0.01, K y = K z = 10, dotted line: K x = 0.01, K y = K z = 1000, and dash-dot line: K x 
= 0.01, K y = K z = 50000. In our energy units Rcj^ = 0.1414, whereas hcu y (= hu) z ) = 4.47, 44.72, and 316.2 for K y = 10, 
1000, and 50000, respectively. The number of bosons (N) used is 10000. The inset shows a comparison between C v 's of an 
anisotropic case (dotted lines for K x = 0.01, K y — K z — 50000) and a pure ID case (solid line) with K x = 0.01. 
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FIG. 2. The single particle DOS for free bosons in an anisotropic harmonic trap vs energy E. Top panel (a): K x = 0.01 and 
K y — K z = 50000; triangles give the numerically calculated exact DOS, dashed line is the DOS obtained from the analytical 
expression in Eq. 1. Bottom panel (b): K x = 0.01 and K y = K z = 10; thin solid line represents numerically calculated exact 
DOS, dashed line is the DOS obtained from Eq. 1. In our energy unit %lo x — 0.1414. 
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FIG. 3. The variation of the specific heat (C v ), bose condensate fraction (No/N), and the fractional population in the 
x-direction (Nm /N) with temperature for bosons in combined optical lattice and anisotropic harmonic potentials. Solid line: 
k x — 0.01, k y = k z — 10, dotted line: k x = 0.01, k y — k z = 100, and dash-dot line: k x = 0.01, k y — k z = 1000 (in units of t 
— 1). The number of bosons (N) used is 10000. 
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FIG. 5. Comparison of the specific heat (C v ) of 10000 bosons in the 3D anisotropic harmonic and optical lattice potentials 
with that of bosons in a ID combined harmonic and optical lattice. Dotted line: k x = 0.01, k y = k z — 1000, dash line: k x = 
0.01, k y = k z — 2000, and dash-dot line: k x — 0.01, k y — k z — 4000 (in units of t = 1). The solid line is for a ID lattice with 
k x = 0.01. 
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FIG. 6. The variation of the specific heat (C v ), bose condensate fraction (No/N), and the fractional population in the 
^-direction (Nm/N) with temperature for bosons in combined isotropic harmonic (k x — k y — k z — 0.01) and anisotropic 
optical lattice potentials. Solid line: t y — t z — 1.0, dotted line: t y = t z = 0.01, and dash-dot line: t y = t z = 0.001 (all energy 
parameters are in units of t x = 1.0). The number of bosons (N) used is 10000. 
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